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Abstract 

We study the multicritical behavior arising from the competition of two dis- 
tinct types of ordering characterized by 0(n) symmetries. For this purpose, we 
consider the renormalization-group flow for the most general 0(rti)©0(ri2)- 
symmetric Landau-Ginzburg- Wilson Hamiltonian involving two fields (fi\ and 
cj>2 with n\ and ni components respectively. In particular, we determine in 
which cases, approaching the multicritical point, one may observe the asymp- 
totic enlargement of the symmetry to O(iV) with iV = m + n^. 

By performing a five-loop e-expansion computation we determine the fixed 
points and their stability. It turns out that for N = n\ + ri2 > 3 the 0(N)- 
symmetric fixed point is unstable. For N = 3, the multicritical behavior 
is described by the biconal fixed point with critical exponents that are very 
close to the Heisenberg ones. For N > 4 and any n\,rii the critical behavior 
is controlled by the tetracritical decoupled fixed point. 

We discuss the relevance of these results for some physically interesting 
systems, in particular for anisotropic antiferromagnets in the presence of a 
magnetic field and for high-T c superconductors. Concerning the SO (5) theory 
of superconductivity, we show that the bicritical 0(5) fixed point is unstable 
with a significant crossover exponent, (^4.4 ~ 0.15; this implies that the 0(5) 
symmetry is not effectively realized at the point where the antiferromagnetic 
and superconducting transition lines meet. The multicritical behavior is either 
governed by the tetracritical decoupled fixed point or is of first-order type if 
the system is outside its attraction domain. 

PACS Numbers: 64.60.Kw, 05.70.Jk, 74.25.Dw, 75.50.Ee 
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I. INTRODUCTION. 



The competition of distinct types of ordering gives rise to multicritical behavior. More 
specifically, a multicritical point (MCP) is observed at the intersection of two critical lines 
characterized by different order parameters. MCP's arise in several physical contexts. The 
phase diagram of anisotropic antiferromagnets in a uniform magnetic field H» parallel to the 
anisotropy axis presents two critical lines in the temperature- H» plane, belonging to the XY 
and Ising universality classes, that meet at a MCP [1],@,[3] • A MCP is also observed in 4 He. It 
arises from the competition of crystalline and superfluid ordering in the temperature-pressure 
phase diagram [|J. MCP's are also expected in the temperature-doping phase diagram of 
high-T c superconductors. Within the SO (5) theory of high-T c superconductivity, it has 
been speculated that the antiferromagnetic and superconducting transition lines meet at a 
MCP in the temperature-doping phase diagram, which is bicritical and shows an effective 
enlarged 0(5) symmetry. On the other hand, the recent experimental evidence of a coex- 
istence region between the antiferromagnetic and superconducting phases is suggestive of a 
tetracritical behavior |7]]. A MCP should also appear in the temperature baryon-chemical- 
potential phase diagram of hadronic matter, within the strong-interaction theory with two 
massless quarks 

Different phase diagrams have been observed close to a MCP. If the transition at the MCP 
is continuous, one may observe either a bicritical or a tetracritical behavior. A bicritical 
behavior is characterized by the presence of a first-order line that starts at the MCP and 
separates the two different ordered low-temperature phases, see Fig. |I[ In the tetracritical 
case, there exists a mixed low-temperature phase in which both types of ordering coexist and 
which is bounded by two critical lines meeting at the MCP, see Fig. 0. It is also possible that 
the transition at the MCP is of first order. A possible phase diagram is sketched in Fig. || 
In this case the two first-order lines, which start at the MCP and separate the disordered 
phase from the ordered phases, end in tricritical points and then continue as critical lines. 

If the order parameters have respectively n\ and n 2 components and the interactions 
are invariant under O(rii) and 0(712), the critical behavior at the MCP can be studied 
by starting from the most general Landau- Ginzburg- Wilson (LGW) Hamiltonian that is 
symmetric under 0(ni)(BO(n 2 ) transformations and contains up to quartic terms 0: 



H 



d d X {^[(d^) 2 + (0^2 



+ 



ri<Pl + r 2 4) + i [u^lf + uMlf + 2w<ft$] } . (1.1) 



Here, the two fields 0i and <f> 2 have n\ and n 2 components respectively. The critical behavior 
at the MCP is determined by the stable fixed point (FP) of the renormalization-group (RG) 
flow when both r% and r 2 are tuned to their critical value. An interesting possibility is 
that the stable FP has 0(N) symmetry, iV = n% + n 2 , so that the symmetry gets effectively 
enlarged when approaching the MCP. This picture has been put forward for the multicritical 
behavior of anisotropic antiferromagnets in an external magnetic field [0,0], for systems with 
quadratic and cubic anisotropy |p!0| , [lT| , [12| , and for high-T c superconductors |5|,|6|, p!3|JT4 1 . 



The phase diagram of the model with Hamiltonian ( |1 . 1|) has been investigated within the 
mean-field approximation in Ref. || (see also Ref. [HJ) 
of a bicritical or tetracritical point, as observed experimentally. 



This analysis predicts the existence 
The nature of the MCP 



depends on the sign of the quantity A = U\U 2 — w 2 , which is relevant in the study of the 
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FIG. 1. Phase diagram in the plane T-g presenting a bicritical point. Here, T is the temperature 
and g a second relevant parameter. The thick line ( "flop line" ) represents a first-order transition. 



stability domain of the Hamiltonian ( |1 . 1| ) . If A > the MCP is tetracritical as in Fig. |2|, 
while for A < it is bicritical, as in Fig. [l]. 

The critical behavior of the model has been investigated in the framework of the e 
expansion 0|§. A low-order calculation HQ shows that the isotropic 0(iV)-symmetric FP 
(JV = rii + 712) is stable for iV < N c = 4 — 2e + 0(e 2 ). With increasing N, a new FP named 
biconal FP (BFP), which has only 0(ni)©0(n2) symmetry, becomes stable. Finally, for 
large N, the decoupled FP (DFP) is the stable FP. In this case, the two order parameters 
are effectively uncoupled at the MCP. The extension of these 0(e) results to three dimensions 
suggests that for n\ — 1 and ri2 = 2, the case relevant for anisotropic antiferromagnets, the 
MCP belongs to the 0(3) universality class, while for n\ = 2 and n<i = 3, of relevance 
for the SO (5) theory of high-T c superconductivity, the stable FP is the BFP. The 0(e) 
computations provide useful indications on the RG flow in three dimensions, but a controlled 
extrapolation to e = 1 requires much longer series and an accurate resummation exploiting 
their Borel summability. As we shall see, the above-reported hypotheses on the three- 
dimensional systems with n% = 1, n-i = 2 and n\ = 2, ri2 = 3 will be both contradicted by a 
higher-order analysis. 

The stability properties of the DFP can be established using nonperturbative arguments 
[ TB| , P~T|JT~8| , P~5| ] , which allow us to compute the RG dimension y w of the operator w<^)\ at the 
DFP. The stability of the DFP depends on the sign of y w : if y w < 0, the DFP is stable. It 
turns out that in three dimensions y w > for N < 3, and y w < for N > 4 for any n\ 
and ri2, showing that the DFP is stable for iV > 4. We should note that the stability of the 
DFP does not allow us to exclude the existence of another stable FP. This possibility, which 
is usually considered rather unlikely flI8| , has been put forward j2(J to explain the Monte 
Carlo results of Refs. |l4|,p0| , which apparently support the stability of a multicritical 0(5) 



FP. 

The phase diagram of the model (|1.1| ) was studied in Refs. 0,0,^]. The DFP is ex- 
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FIG. 2. Phase diagram with a tetracritical point. 



pected to be generically tetracritical: indeed, in this case the MCP should correspond to 
a generic intersection of the two critical lines with O(ni) and 0(712) symmetry. The sta- 
ble O(N) FP — as we shall see, this is the case only for N = 2 — can be either bicritical or 
tetracritical. The possibility of two different phase diagrams for the same FP is due to the 
presence of a dangerously irrelevant operator . Little is known for the BFP, although 

a phenomenological extension of the mean-field arguments would predict a tetracritical be- 
havior [^] . When the initial parameters of the Hamiltonian are not in the attraction domain 
of the stable FP, the transition between the disordered and ordered phases should be of 
first order in the neighborhood of the MCP [EME3LE3J. However, the transition along the 



critical lines may become continuous sufficiently far from the MCP |f25| , p26| . A possible phase 
diagram is sketched in Fig. |3|. 

In this paper we extend the analysis of the multicritical RG flow to 0(e 5 ). The stability 
of the O(N) FP is also discussed in the framework of fixed- dimension expansion in three 
dimensions, for which six-loop series have been computed. These calculations allow us to 
obtain a rather conclusive picture of the multicritical RG flow in three-dimensional systems. 
In particular, the O(N) FP is stable only for N = 2. Therefore, the symmetry enlargement 
occurs only when the competing order parameters have Ising symmetry. For N > 3, the 
O(N) FP is unstable and therefore the enlargement of the symmetry to O(N) at the MCP 
requires an additional tuning of the parameters: beside tuning r\ and r 2 , a third parameter 
must be properly fixed to decouple the additional relevant interaction. The crossover expo- 
nent associated with this RG instability is 04,4 ~ 0.01 for iV = 3, ^4,4 ~ 0.08 for iV = 4, 
04,4 ~ 0.15 for iV = 5, and 4 ,4 1 for iV — > 00. For iV = 3 the stable FP is the BFP. The 
critical exponents are however very close to the Heisenberg ones, so that distinguishing ex- 
perimentally the 0(3) FP and the BFP is a very hard task, taking also into account the very 
small crossover exponent governing the unstable flow from the 0(3) FP. The case N = 5, 
n\ = 2, n 2 = 3 is relevant for the SO (5) theory |5]|| of high-T c superconductors, which 
proposes a description in terms of a three-component antiferromagnetic order parameter 
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FIG. 3. Phase diagram with a first-order MCP. The thick lines represent first-order transitions. 



and a c?-wave superconducting order parameter with U(l) symmetry, with an approximate 
0(5) symmetry. For N = 5 the only stable FP is the DFP which predicts, if the transition 
is continuous, a tetracritical behavior. This may explain a number of recent experiments, 
see, e.g., Refs. |]7|,^, PB| , p9| , |50] , |3"T] ,|5^| , that provided evidence of a coexistence region of the 
antiferromagnetic and superconducting phases. The 0(5) FP is unstable with a crossover 
exponent ^44 ~ 0.15, which, although rather small, is nonetheless sufficiently large not to 
exclude the possibility of observing the RG flow towards the eventual asymptotic behavior 
for reasonable values of the reduced temperature even in systems with a moderately 
small breaking of the 0(5) symmetry, for instance in those described by the projected SO (5) 
model discussed in Refs. |||34,35]. Of course, when the effective Hamiltonian parameters 
are outside the attraction domain of the stable FP, the transition at the MCP is expected 
to be of first-order type. Some of the results concerning the stability properties of the 0(N) 



FP were already presented in Ref. . 

The paper is organized as follows. In Sec. |I| we present our five-loop calculations in the 
framework of the e expansion. In Sec. |TTI|we discuss the stability of the 0(iV)-symmetric FP 
under generic perturbations. The results are then applied to establish the stability properties 
of the O(N) FP. In Sec. f[V] the multicritical RG flow is analyzed. In Sec. [V] we draw our 
conclusions and discuss their relevance for some physical systems. 



II. e EXPANSION OF THE 0(m)ffiO(n 2 ) THEORY 

We extended the e expansion of the critical exponents at the different FP's for the 
0(ni)©0(n2) symmetric theory to 0(e 5 ). For this purpose, we considered the minimal 
subtraction (MS) renormalization scheme [37]. We computed the divergent part of the 
irreducible two-point functions of the fields </>i and 02, of the two-point correlation functions 
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with insertions of the quadratic operators and 2 , and of the three independent four- 
point correlation functions (0i • 0i 0i ■ 0i), (0i • 0i 2 • 2 ), and (0 2 '02 02 • 02)- The 
diagrams contributing to this calculation are a few hundreds. We handled them with a 
symbolic manipulation program, which generated the diagrams and computed the symmetry 
and group factors of each of them. We used the results of Ref. P5[ , where the primitive 
divergent parts of all integrals appearing in our computation are reported. We determined 
the renormalization constants and associated with the fields 0i and 2 respectively, 
the 3x3 renormalization matrix Zfj of the quartic couplings defined by gs,i = fi e ZfjgRj 

where gs,i = («i, w 2 , w), and the 2x2 renormalization matrix Zfj of the quadratic operators 

<p\ and 2 . The (3 functions Pi(gRj) and the RG dimensions 7^, 7^3, 7^ are determined 
using the relations 

dgRA 



9/i 



(2.1) 



9B,j 

7*0aw) = E#a^, (2 ' 2) 

7gW = £/^(^V" (2-3) 

The zeroes g}^ j of the /9-functions provide the FP's of the theory. In the framework of the 
e expansion, they are obtained as perturbative expansions in e and are then inserted in the 
RG functions to determine the e expansion of the critical exponents. The stability of each 
FP is controlled by the 3x3 matrix 



_ dPiig^k) 
d 9RJ 



(2.4) 

9R,k=g* Rk 



The two exponents 7/1 and 7/ 2 , related to the short-distance behavior of the two-point func- 
tions of the fields 0i an d 02, are given by 771 = 70 1 (g|jj) and 7/ 2 = l^ 2 i.9Ri)- From the 
eigenvalues v\ and z/ 2 of the matrix 7^ , if V\ > v% y one obtains v = v\ and = ^i/V 2 , where 
is the crossover exponent associated with the quadratic instability. 

We performed several checks of the perturbative series. In particular, the critical- 
exponent series agree with the existing 0(e 5 ) ones for the 0(iV)-symmetric theory |3T||4"(1 



in 



the proper limit. Morover, as we shall discuss in the following section, we can also compare 
with some results for the O(N) theory in the presence of cubic anisotropy [IT|, finding agree- 



ment. Some of the five- loop perturbative series will be reported in the following sections. 
The complete list of series is available on request. 

Since the e expansion is asymptotic, the series must be properly resummed to provide 
results for three-dimensional systems. We used the Pade-Borel method except for the series 
at the O(N) FP. In this case, we applied the conformal- mapping method that takes into 



account the known large-order behavior of the expansion. See, e.g., Refs. [[43[[44[ for reviews 
of resummation methods. 
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III. STABILITY OF THE O(N) FIXED POINT 



In this section we discuss the stability of the 0(A) FP, where A = n\ + n 2 , to establish 
in which cases the enlargement of the symmetry 0(ni)©0(ri2) to 0(A) is realized at the 
MCP without the need of further tunings. 

Let us consider the general problem of an 0(A)-symmetric Hamiltonian in the presence 
of a perturbation P, i.e., 

H = j d d x [\(d^f + ±r$ 2 + ±u($ 2 ) 2 + h p P] , (3.1) 

where $ is an iV-component field and h p an external field coupled to P. Assuming P to be 
an eigenoperator of the RG transformations, the singular part of the Gibbs free energy for 
the reduced temperature t — > and h p — > can be written as 

^sing(*,M« \t\ dv ?(h p \t\-^), (3.2) 

where 4> p = y p v is the crossover exponent associated with the perturbation P, y p is the RG 
dimension of P, and T{x) is a scaling function. If y p > the pertubation is relevant and its 
presence causes a crossover to another critical behavior or to a first-order transition. 

In order to discuss the stability of the 0(A) FP in general, we must consider any per- 
turbation of the 0(A) FP. We shall first consider perturbations that are polynomials of the 
field $ a . Any such perturbation can be written |^5| as a sum of terms P^l'"' at , m > I, which 



are homogeneous in $ a of degree m and transform as the Z-spin representation of the 0(A) 
group. Explicitly, we have 

p«i.~.«« = (<$>2}m-lQai,:, ai (3.3) 

where g^'-« a « is a homogeneous polynomial of degree I that is symmetric and traceless in 
the I indices. The lowest-order even polynomials are 

Qab = $ a $ 6 _ }_ 6 ab^2 (3 ^ 

Qabcd _ ^a^b^c^d _ _L_<|) 2 (jj ab § c § d _|_ $ ac $ b $ d _|_ fi ad $ b $ c _|_ ^ 6c $ a $ d _|_ ^ bd $ a $ c _|_ ,5 cd $ a $ b 
(V+2)(V+4) 



+ 77vTWlvT4T( $2 ) 2 ( SabSCd + + 5a ^" C • ( 3 - 5 ) 



The classification in terms of spin values is particularly convenient, since polynomials with 
different spin do not mix under RG transformations. On the other hand, operators with 
different m but with the same / do mix under renormalization. At least near four dimensions, 
we can use standard power counting to verify that the perturbation with indices m, / mixes 
only with P m ',ii w! < m. In particular, Pij renormalizes multiplicatively and is therefore a 
RG eigenoperator. Moreover, if y mt i is the RG dimension of the appropriately subtracted 
P m ,u one can verify that for small e, y m j < 0, for I > 5, i.e. the only relevant operators 
have I < 4. We will assume this property to hold up to e = 1. We notice that it is 
certainly incorrect in two dimensions where perturbations are relevant (A > 3) or marginal 



(A = 2) for all values of /. j|6| In principle, we should also consider terms with derivatives 



of the field, but again, using power counting, one can show that they are all irrelevant or 
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redundant. Therefore, beside the 0(iV)-symmetric terms $ 2 and ($ 2 ) 2 there are only three 
other perturbations that must be considered, P$ b 2l P£ b 2 , and Pf A d . Note that, according 
to the above-reported discussion, P 2 \ and Pff d are RG eigenoperators, while P%\ must 
be in general properly subtracted, i.e. the RG eigenoperator is PI \ + zP 2 \ for a suitable 
value of z. The determination of the mixing coefficient z represents a subtle point in the 



fixed-dimension expansion |47[], but is trivial in the MS scheme in 4 — e dimensions, in which 
operators with different dimensions never mix so that z = 0. 

According to the above-presented general analysis, the stability properties of the O(N) 
FP can be obtained by determining the RG dimensions of the five operators reported above. 
Of course, the result does not depend on the specific values of the indices and thus one 
can consider any particular combination. We now show that such dimensions determine the 
crossover exponent <fi and the eigenvalues of the stability matrix Q at the O(N) FP for the 
0(ni)©0(ri2) theory. Starting from the general expressions, one can construct combinations 
that are invariant under the symmetry group 0(ni)©0(n 2 ). Explicitly, they are given by 

a=l iV 
P 4 ,0 = ($ 2 ) 2 , 7>4,2 = &V 2 ,2, 

^4,4 = 2^ 1^ P 4,4 = 0102 + (AT , oV » r , A y (3-6) 



o=l b=ni+l 



iV + 4 (iV + 2)(iV + 4)' 



Here $ is the iV-component field (0i,02)- The RG dimensions of P2,o an d of V^o are well- 
known and can be computed directly in the 0(iV)-invariant theory. In particular, y 2 ,o = 1/^ 
and 2/4^ = — to, where uo is the leading irrelevant exponent in the O(N) -invariant theory. 
The RG dimension 7/2,2 of 7^2,2, and therefore of the operator P^, provides the crossover 
exponent (f> = yipV at the MCP. We denote such exponent by 4>t to stress the fact that it is 
associated with the tensor quadratic operator. Setting 



1 + J>e\ (3.7) 



we obtain at five loops 



N _ Af(Af 2 +24Af+68) 

2(N+8) ' t' 2 ~ 4(Af+8) 3 



P3 



iV(iV 4 +48iV 3 +788iV 2 +3472iV+5024) _ 6Af(5Af+22)C(3) 
8(Af+8) 5 (Af+8) 4 ' 



_ Af(Af 6 +72Af 5 +2085A r4 +28412Af 3 +147108Af 2 +337152Af +306240) . iy(-iy 4 +13jV 3 -544iV 2 -4716iV-8360K(3) 
^ 4 ~~ 16(Af+8) 7 + (Af+8) 6 



P5 



N(5N+22)n 4 . 20Af(2Af 2 +55Af+186)C(5) 
20(7V+8) 4 r (Af+8) 5 

N (17677824+28388096 jV+19390624 jV 2 +6723904 jV 3 +1177480 Af 4 +95668 jV 5 +4154 7V 6 +96 N 7 +N 8 ) 

32 (8+AT) 9 

N (8360+4716 7V+544 Af 2 -13 N 3 +N 4 ) tt 4 . 5 N (186+55 N+2 N 2 ) tt 6 

120(8+A r ) 6 ' 189(8+A r ) 5 

N (554064+465592 Af+125232 jV 2 +7584 jV 3 -661 jV 4 +9 Af 5 ) ((3) 
(8+Af) 8 

2 N (24528+14468 7V+2028 Af 2 +39 Af 3 +4 N 4 ) C(3) 2 
(8+Af) 7 

N (466016+280596 Af+33832 Af 2 -2857 Af 3 -230 Af 4 ) <(5) _ 441 N (526+189 N+U N 2 ) ((7) /o o\ 

2 (8+Af) 7 2(8+Af)6 • 
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TABLE I. Estimates of the RG dimensions y 2 2, 2/4,0) 2/4, 2, and 2/4,4, and of the crossover 
exponents 4>t = 2/2,2^, 04,4 = 2/4,4^ as obtained by various approaches: e expansion (e exp), 
fixed-dimension expansion (d = 3 exp), high-temperature expansion (HT exp), Monte Carlo simu- 
lations (MC), and 1/N expansion (1/N exp). Their values in the large-iV limit, see, e.g., Ref. [51], 
are also reported. 
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n 077/ 
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El 
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u.uuy^4 j 




ITT OVT-> 

n 1 exp 




1 9^0/1 IK"' 
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-0.773 1 58] 






-0.0007(29) |3 






1/N exp 
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e exp 


1.813(6) 


1.329(16) 


-0.795(30) 


y 


-0.493(14) 


0.105(6) || 


0.079(5) 




d = 3 exp 




1.35(4) @ 


-0.774(20) 






0.111(4) H 


0.083(3) 




MC 






-0.765 |H 






0.130(24) |3 






1/iV exp 
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5 


e exp 


1.832(8) 


1.40(3) 


-0.783(26) 




-0.441(13) 


0.198(11) 


0.151(9) 




d = 3 exp 




1.40(4) H 


-0.790(15) 






0.189(10) H 


0.144(8) 




MC 




1.387(30) {D 


1 












1/JV exp 




1.422 §§ 
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This series extends the three-loop results of Ref. |48| and the four-loop results of Ref. [^9 
In the appropriate limit, it is in agreement with the 0(N~ 2 ) expression of Ref. ||50|| . In 
Table [I] we report the estimates of 1/2,2 and 4>t for N = 2, 3, 4, 5 obtained from the analysis 
of the five-loop perturbative expansion (|3.8|). As expected, since 2/2,2 > in all cases, the 
quadratic perturbation P 2 a 2 I s always relevant. The results are compared with the estimates 
obtained from the analysis of its six-loop fixed- dimension expansion 13 and of its large- A" 



expansion to 0(1/N 2 ) ||50|| , and by using high-temperature techniques []55] and Monte Carlo 
simulations [|14[]. We also mention that consistent results were obtained from the analysis 
of the four-loop series of < 
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1.177 for N = 2 and 



_ = 1.259 for N = 3. Some 
experimental results for 0V can be found in Ref. [4~3 . 

The perturbative expansions of the RG dimensions of the operators V^i, and therefore of 
the more general operators P^i, can be obtained from the eigenvalues of the stability matrix 
Q at the O(N) FP. For this purpose, it is convenient to perform a change of variables, 
replacing u\, «2, and w with gi, I = 0,2,4, which are the quartic couplings associated with 
the operators V m 1 and are explicitly defined by the relation 



9 2 f + 



90^,0 + #2^4,2 + 04^*4,4 ■ 



(3.9) 



In this basis Q is diagonal and the eigenvalues of Q are simply given by 
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d{3i{g , g 2 , g*) 

LOi = - 



dgi 



(3.10) 

90=9* N , 92=0, 34=0 



where / = 0,2,4, fii are the /3-functions associated with the couplings gi, and g* N is the FP 
value of the quartic coupling in the 0(A^)-symmetric theory. The critical exponent uj q is the 
leading irrelevant operator in the 0(iV)-symmetric theory. Its 0(e 5 ) expansion can be found 
in Refs. [|39| ) fl0|l ; several estimates are reported in Refs. [|52],[44|]. The RG dimension y^i of 
the perturbation P^i is given by 

y i} i = -oj h (3.11) 
We report here the five-loop e expansion of and 2/4 4. Setting 

2/4,/ = E c Me\ (3-12) 

i=l 

we have 

„ _ 8 „ _ 336+68 N+7N 2 

c 2,l - N+s , (-2,2 - ( 8 +AT)3 , 

76544+26176 jV+3264 Af 2 +28 N 3 -N 4 _ 12 (352+82 A^+7Af 2 ) C(3) 
C 2,3 4(8+Af) 5 (8+Af) 4 

_ (20796416+10251520 iV+2207744 Af 2 +271328 Af 3 +24824 jV 4 +820 iV 5 +5 N 6 ) _ (352+82 N+7 N 2 ) tt 4 
C 2,4 — 16(8+iV) 7 10(8+Af) 4 

2 (-92928-34776 JV-7544 AT 2 -928 Ar 3 -67 N 4 +N 5 ) C(3) , 80 (2232+632 iV+60 N 2 +N 3 ) C(5) 
(8+N) 6 (8+iV) 5 

_ 6019366912+3720851456 ^+994704384 ^+135243264 iV 3 +6891584 Af 4 -590816 Af 5 -60520 V 6 -1732 Af 7 -13 N 8 
C 2,5 — 64 (8+ AO 9 



, (92928+34776 7V+7544 Af 2 +928 A r3 +67 N 4 —N 5 ) n 4 , 20 (2232+632 jV+60 Af 2 +jV 3 ) tt 6 
' 60(8+V) 6 ' 189(8+Af) 5 

(117872640+62925184 V+14334912 Af 2 +1577392 jV 3 +67848 Af 4 -1872 Af 5 +200 N 6 -7 N 7 ) C(3) 

4(8+Af) 8 

, 8 (104832+100312 iV+24994 AT 2 +2571 Af 3 +83 Af 4 +2 N 5 ) C(3) 2 
+ {8+N) 7 

(7263744+3733728 Af+1095516 jV 2 +170284 jV 3 +14035 A^ 4 +322 N 5 ) C(5) _ 441 (16832+5590 jV+631 jV 2 +23 A f3 )C(7) (o 

(8+N) 7 (8+N) 6 y 6 

and 

_ AT-4 „ _ 152+14 N+5N 2 

C 4,l — TTTaF; C 4,2 



8+JV' ^> z (8+N) 3 ' 

17024-1568 jV-1464 Af 2 -398 N 3 -13N 4 _ 48 (46+7 N+N 2 ) C(3) 
4(8+Af) 5 (8+iV) 4 ' 

(2995712+402304 Af+223328 jV 2 +112856 iV 3 +27272 Af 4 +1516 jV 5 +29 A" 6 ) 2 (46+7 N+N 2 ) tt 4 



C4 > 4 16(8+Af) 7 5(8+iV) 



c 4,5 



3 (-21568+1664 V+1592 jV 2 +256 Af 3 -8 Af 4 +Af 5 ) C(3) , 120 (712+130 AT+13 N 2 ) <(5) 
(8+Af) 6 (8+Af) 5 ' 

-365813760+95377408 A r +75546624 A r2 +35042816 A r3 +11477472 Af 4 +2184488 A r5 +148600 AT 6 +4712 A r7 +61 N s 



64(8+AT) 
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(-21568+1664 Af+1592 Af 2 +256 N 3 -8 N 4 +N 5 ) n 4 , 10 (712+130 AT+13 N 2 ) tt 6 
40 (8+ AT) 6 ' 63 (8+AO 5 

(37827072+13773568 N +3633344 A^ 2 +689728 Af 3 +54184 Af 4 -3272 Ar 5 +188 Af 6 -5 V 7 ) C(3) 

4 (8+Af) 8 

, 12 (11456+24112 Ar+6648 A r2 +790 Af 3 +5 N 4 ) C(3) 2 
+ (8+AT) 7 

2 (1018944+128152 AT-3060 A 1 " 2 — 9018 A" 3 +347 Af 4 +12 A" 5 ) ^(5) _ 2646 (1268+272 AT+25 N 2 +N 3 ) ((7) /o i a\ 

(8+N) 7 (8+N) 6 ■ V 6 - iq: ) 
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At one loop, these results agree with those reported in Ref. | 45]| . The results of the analyses 
of these series are reported in Table [p. They show that y 4i2 is always negative, so that the 
corresponding spin-2 perturbation Pf^ is always irrelevant. On the other hand, the sign of 
1/4,4 depends on iV: it is clearly negative for N = 2 and positive for N > 4. For N = 3 it is 
marginally positive, suggesting the instability of the 0(3) FP. This fact will be confirmed by 
the more accurate results discussed below. The corresponding crossover exponents m) ; = 
y m} i v can be determined using the following estimates of v: v = 0.67155(27) for N = 2 [55 
v = 0.7112(5) for N 



3 59 



v 



0.749(2) for N = 4 f58[, and v = 0.762(7) for iV = 5 ]60 



Other estimates of v can be found in Ref. 

The RG dimension 2/44 can also be obtained starting from the cubic-symmetric LGW 
Hamiltonian, 



r f 1 N 1 N N ) 

J d d x ~ E [( + + ~ [«(£ + « £ «>fl , (3-15) 

I i=l i i ) 



see, e.g., Ref. [16], and in particular from the results for the stability properties of the O(N) 
FP in the presence of a cubic-symmetric anisotropy. The point is that the cubic-symmetric 
perturbation is a particular combination of the spin-4 operators P£ b £ d and of the spin-0 term 
($ 2 ) 2 . Indeed, one may rewrite 

N N 

E(f)4 _ \ ^ paaaa , 
{ 2^ 4 ' 4 jV + 2 

i=l a=l iV -r ^ 



($ 



,2\2 



(3.16) 



Thus, the stability of the 0(N) FP against the cubic-symmetric perturbation J2i&t 1S con- 
trolled by the RG dimension 2/4 4 of the spin-4 operator P^ d . 

The RG flow for the cubic-symmetric theory has been investigated by employing field- 
theoretical methods, based on perturbative expansions |^ , |5^ , |6l] , |6^ , |6^ , |6^ , |65| , |66|[ or approxi- 
mate solutions of continuous RG equations |67|J6^|J69]] , and lattice techniques, such as Monte 
Carlo simulations |57] and high-temperature expansions JF(J; see, e.g., Ref. [Q for a recent 
review. In particular, the RG functions have been computed to five loops in the e expan- 
sion |41]] , and to six loops in a fixed- dimension expansion in powers of the zero- momentum 
quartic couplings [|k|. In these perturbative schemes 



2/4,4 



dp v (u,v) 



dv 



(3.17) 



u=g* N , v=0 



where (3 V is the /3-function associated with the quartic coupling v, and g* N is the FP value of 
the quartic coupling in the 0(iV)-symmetric theory. This allows us to determine 7/^4 using 
the five-loop expansions reported in Ref. ^T|. We reobtain again Eq. ( |3.14j ), confirming 
the correctness of our calculation. Moreover, using Eq. (|3.17|) and the results reported 
in Ref. |p3|| , one can also compute 7/44 in the framework of the fixed- dimension expansion 
to six loops. The resulting estimates, obtained by using the conformal-mapping method, 
are reported in Table |. They show that the spin-4 perturbation P 4 4 is relevant for all 
iV > 3. In Table | the Monte Carlo results of Ref. [57] are also shown; they were obtained 
by simulating the standard iV-vector model and computing the RG dimension of the cubic- 
symmetric term Ya s ti where Si is the A^-component spin variable. We may also consider 
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the value N c such that for N > N c the cubic-symmetric anisotropy, and therefore the spin-4 
perturbation P 44 , becomes relevant at the 0(A) FP. All studies reported in the literature 
indicate A c « 3 and definitely N c < 4, see, e.g., Refs. |5|j43J6gggj6|j63J6|j6|j53J7C 



The most accurate results have been provided by analyses of high-order perturbative field- 
theory expansions, which predict iV c < 2.9 in three dimensions. In particular, different 
analyses of the six- loop fixed- dimension series yielded the estimates A c = 2.89(4) |53| and 
N c = 2.862(5) [p^] ; similar results were also obtained from shorter series, see, e.g., Refs. 
|63| , |64]| . These results have been confirmed by the analysis of the 0(e 5 ) series []4"1], [66| , |5"3|1 . A 



constrained analysis taking into account the two-dimensional value of A c , N c = 2, provided 
the estimate N c = 2.87(5) |53|], which makes the evidence supporting y 4)4 > for A = 3 
stronger than the estimate j/ 4)4 = 0.003(4) obtained from the direct analysis of its 0(e 5 ) 
series. 

In conclusion, these results provide a rather robust evidence that for A > 3 the 0(A) 
FP is unstable with respect to spin-4 perturbations P 4 4 cd , and, as a consequence, that the 
0(A) FP is a unstable MCP for iV > 3. 



IV. RG FLOW AT THE MULTICRITICAL POINT 

As already shown by the 0(e) computations of Ref. ||, the O(ni) ©0(772) theory at the 
MCP has six FP's. Three of them, i.e. the Gaussian, the O(ni) and the 0(77,2) FP's, are 
always unstable. The other three FP's are the 0(A) fixed point (FP), the biconal fixed 
point (BFP), and the decoupled fixed point (DFP). The stability of these FP's depends on 
n 1 and 772. In particular, in the preceding section we have established that the O(A) FP is 
stable for A = 2 and unstable for A > 3, for any n% and 772. 

The stability properties of the DFP can be determined using nonperturbative scaling 
arguments [[T^,[n],[n^[nj . At the DFP, the quartic coupling term wcpfcpl scales as the product 
of two energy-like operators, which have RG dimensions (1 — a^jvi where a, and i>\ are the 
critical exponents of the 0(t7j) universality classes. Therefore, the RG dimension related to 
the w-perturbation is given by 

yu> = -z- + ^- = - + --d. (4.1) 

Note that this relation is satisfied order by order in the e expansion. Indeed, the e expansion 
of y w obtained from the stability matrix Q at the DFP coincides with the series obtained 
from the right-hand side of Eq. ( ft.lfl , using the five-loop expansions of for the O(rij) 
universality classes. Taking into account that the DFP is stable with respect to the other 
two RG directions, one can determine the stability properties of the DFP from the sign of 
y w . Using the estimates of the critical exponents of the three-dimensional O(rij) universality 



classes (see, e.g., Ref. [§4| for a review), y w turns out to be negative for iV = ri\ + 772 > 4, 
and positive for iV = 2, 3. [I7T[] Three-dimensional estimates of y w for N < 5 are reported in 
Table [Q]. These results show that the tetracritical DFP is stable for N > 4 for any n 1; ra 2 . 

The results concerning the O(A) FP and the DFP suggest that the stable FP for N = 3 
is the BFP. This is substantially confirmed by the five-loop analysis of the stability matrix Q 
at the BFP. Below we report the expansions of the critical exponents at the BFP for n\ = \ 
and 772 — 2. The eigenvalues of the stability matrix Q are 
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TABLE II. Estimates of the RG dimension y w at the DFP. They are obtained by using Eq. ( |4.l| 
and the estimates of the O(iV) critical exponent u reported in Refs. [72,75,56,59|,58|. 

N = n-y + n 2 



2 
3 
4 



~T 
l 
l 

2 
1 

2 



"2 

"T 

2 
3 
2 
4 
3 



0.1740(8) 

0.0761(7) 

-0.0069(11) 

-0.0218(12) 

-0.078(4) 

-0.1048(12) 



WMjl = e - 0.579364 e 2 + 1.344815 e 3 - 4.058162 e 4 + 14.526420 e 5 + 0(e 6 ), 

w b i,2 = 0.491105 e - 0.084149 e 2 + 0.361174 e 3 - 0.776741 e 4 + 2.593212 e 5 + 0(e 6 ), 

Wbi3 = -0.130195 e + 0.278782 e 2 - 0.379711 e 3 + 0.868886 e 4 - 2.656984 e 5 + 0(e 6 ). (4.2) 

The expansions of the critical exponents are 

?7 bi(1 = 0.0208306 e 2 + 0.0182941 e 3 - 0.00777325 e 4 + 0.0210296 e 5 + 0(e 6 ), 

77bi, 2 = 0.0200806 e 2 + 0.0195184 e 3 - 0.00848020 e 3 + 0.0236398 e 5 + 0(e 6 ), 

v^ x = \ + 0.1114875 e + 0.0667684 e 2 - 0.00616190 e 3 + 0.0779498 e 4 - 0.193367e 5 + 0(e 6 ), 

^bi, 2 = \ + 0.0234143 e + 0.0289670 e 2 - 0.00547548 e 3 + 0.0381483 e 4 - 0.106076 e 5 + 0(e 6 ), 

bi = ^iii = 1 + 0.176147 e + 0.0673541 e 2 - 0.0147318 e 3 

UU ' 2 + 0.0783198e 4 - 0.190099 e 5 + 0(e 6 ). (4.3) 

We analyzed these series using the Pade-Borel resummation method. The estimates of the 
eigenvalues of the stability matrix are u^i = 0.79(2), u>u,2 = 0.57(4), and co> bij3 = 0.01(1). 
They are all positive, supporting the stability of the BFP, although the result for co> bi 3 is not 
sufficiently precise to definitely exclude the opposite sign. Concerning the critical exponents, 
we obtained r/u,i = 0.037(5), r/u,2 = 0.037(5), = z/bi,! = 0.70(3), and (fiu = 1.25(1). Note 
first that 77^,1 ~ Vbi,2, as it can be directly guessed by looking at the coefficients of their 
expansions. A direct analysis of their difference gives the bound \rjuA — VbiM ~ 0.0005. 
Second, note that, within the errors, the BFP exponents are very close to the Heisenberg 
ones, whose best estimates are t]n = 0.0375(5), ujj = 0.7112(5), and 4>h = 1.250(15) from 
high-temperature techniques p9| , |55| , r\n = 0.0355(25), uh = 0.7073(35), and 4>h = 1.27(2) 



from the six-loop fixed- dimension expansion [p2| , p4[ , t)h = 0.0375(45), vh = 0.7045(55), 
and <pH — 1.260(11) from the five-loop e expansion Rather stringent bounds on the 



differences between the biconal and Heisenberg exponents can be obtained by considering 
the expansions of their differences, which have much smaller coefficients. Their analysis 
yields 

- Vh\ Z 0.0005, 
|v?bi,2 - < 0.0001, 
0.001, 

\<l>u-M < 0.005. (4.4) 
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We have also studied the stability of the BFP for larger values of N. For N = 4, and in 
both cases n = 1, n 2 = 3 and n\ = n 2 = 2, the five- loop calculation gives the expansions 
of the critical exponents at the BFP only to 0(e 4 ), because of the additional degeneracy of 
the 0(4) FP and of the BFP at 0(e). In particular, for the smallest eigenvalue we obtain 

uta, 3 (m =l,na = 3) = |e 2 -0.3306439 e 3 + 0.7376491 e 4 + 0(e 5 ), 

cj bi)3 (m = 2,n 2 = 2) = |e 2 - 0.319872 e 3 + 0.696458 e 4 + 0(e 5 ). (4.5) 

It is difficult to extract reliable estimates from these series. In both cases, we find that u;bi,3 
is small, but we are unable to determine reliably its sign. 

For iV > 5 we find that the BFP is unstable for all values of n\ and n 2 . In particular, 
for N = 5, Tii = 2, n 2 = 3, for the smallest eigenvalue we obtain 

Wbi;3 = 0.052584e + 0.0331401 e 2 - 0.242179 e 3 + 0.358964 e 4 - 1.242100 e 5 + 0(e 6 ), (4.6) 

which gives co>bi,3 = —0.07(5). 



V. CONCLUSIONS AND DISCUSSION 

We have studied the multicritical behavior at a MCP, where two critical lines with O(rzi) 
and 0(n 2 ) symmetry meet. It has been determined by studying the RG flow of the most 
general 0(ni)©0(n2)-symmetric LGW Hamiltonian involving two fields <fri and <p 2 with n\ 
and n 2 components respectively. We have extended the e expansion of the critical exponents 
and of the stability matrix of the FP's, previously known to one-loop order, to five loops. 
The stability of the O(N) FP has also been discussed in the framework of the fixed- dimension 
expansion in three dimensions to six loops. 

The main properties of the RG flow of the 0(n 1 )©0(ra 2 )-symmetric system at the MCP 
can be summarized as follows. 

• The O(N) FP is stable only for N = 2, i.e. when two Ising-like critical lines meet. 
It is unstable in all other cases, i.e. for all ni and n 2 such that n\ + n 2 = N > 3. 
Beside being unstable with respect to the spin-0 and spin-2 quadratic perturbations, 
for N > 3 the O(N) FP is also unstable with respect to quartic perturbations belonging 
to the spin-4 representation of the O(N) group, cf. Eq. ( |3.5| ). This implies that for 
N > 3 the enlargement of the symmetry 0(ni)®0(n 2 ) to O(N) requires an additional 
parameter to be tuned, beside those associated with the quadratic perturbations, r\ 
and r 2 in the LGW Hamiltonian. The associated crossover exponents 04,4 = 7/4,4 v are: 
04,4 « 0.01 for iV = 3, 04,4 « 0.08 for iV = 4, 4 ,4 ~ 0.15 for iV = 5, and 4 ,4 1 for 
N—> 00 (see Table |). 

• For = 3, i.e. for n\ = 1 and n 2 = 2, the critical behavior at the MCP is described by 
the BFP, whose critical exponents turn out to be very close to those of the Heisenberg 
universality class, see Eq. 



For N > A and for any m > 1 and n 2 > 1, the tetracritical DFP is stable. This 
has been inferred using nonperturbative arguments [16. T7. 18, Tj| that allow us to write 
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the relevant stability eigenvalue y w in terms of the critical exponents of the 0(n;) 
universality classes, cf. Eq. (|4.1| ). The e-expansion analysis shows that the BFP is 
unstable for all cases with N > 5, while it is not conclusive for the cases with N = 4. 

• When the initial parameters of the Hamiltonian are not in the attraction domain of 
the stable FP, the transition between the disordered and ordered phases should be of 
first order in the neighborhood of the MCP. In this possible phase diagram is 

given in Fig. [5| Close to the MCP all transition lines are first-order ones. However, far 
from the MCP, the high-temperature transitions may become continuous, belonging 
to the O(ni) and 0(77-2) universality classes. 

As already mentioned in the introduction, a multicritical behavior has been observed in 
several systems. 

Anisotropic antiferromagnets in a uniform magnetic field H\\ parallel to the anisotropy 
axis present a MCP in the T—H\\ phase diagram, where two critical lines belonging to the XY 
and Ising universality classes meet . The results presented above predict a multicritical 
BFP. The mean- field approximation assigns a tetracritical behavior to the MCP 0, but a 
more rigorous characterization, that requires the computation of the corresponding scaling 
free energy, is needed to draw a definite conclusion. Experimentally, the MCP appears to be 
bicritical, see, e.g., the experimental results of Refs. |74] , |75|1 ; numerical Monte Carlo results 
hint at the same behavior, although with much less confidence |[76|| . Our results contradict 
the 0(e) calculations of Refs. suggesting the stability of the 0(3) FP. Notice that it 

is very hard to distinguish the biconal from the 0(3) critical behavior. For instance, the 
correlation-length exponent v differs by less than 0.001 in the two cases. However, one may 
still hope to distinguish the two FP's by measuring some universal amplitude ratio that 
varies more significantly in the two cases. The crossover exponent describing the crossover 
from the 0(3) critical behavior is very small, i.e., 04,4 ~ 0.01, so that systems with a 
small effective breaking of the 0(3) symmetry cross very slowly towards the biconal critical 
behavior or, if the system is outside the attraction domain of the BFP, towards a first-order 
transition; thus, they may show the eventual asymptotic behavior only for very small values 
of the reduced temperature. 

Isotropic antiferromagnets in a magnetic field are quite a special case. Indeed, the critical 
transition at H = is a MCP with 0(3) symmetry, as observed experimentally, see, e.g., 
Ref. [1771]. As we discussed, for H 7^ 0, two relevant perturbations are switched on, and they 
can give rise in principle to a more complex phase diagram. Finally, it should be noted 
that in real antiferromagnets additional nonisotropic interactions are present, giving rise 



to lower-symmetry MCP's. In Ref. |78j the magnetic phase diagram of NiCl2-4H20 was 
studied. The orthorombic symmetry of the crystal gives rise to Ising transition lines both 
for small and large H\\, so that n\ = n% = 1. As predicted by the theory, the MCP is a 
bicritical XY point. A similar experiment is reported in Ref. [f79| . A tetracritical XY MCP 



is observed in anisotropic antiferromagnets when the magnetic field is perpendicular to the 
symmetry axis, see, e.g., Ref. |74[] for an experimental study. 



High-T c superconductors are other interesting physical systems in which MCP's may 
arise from the competition of different order parameters. At low temperatures these ma- 
terials exhibit superconductivity and antiferromagnetism depending on doping. The SO (5) 
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theory |5]|| attempts to provide a unified description of these two phenomena, involving 
a three-component antiferromagnetic order parameter and a d-wave superconducting order 
parameter with U(l) symmetry, with an approximate 0(5) symmetry. This theory predicts 
a MCP arising from the competition of these two order parameters when the corresponding 
critical lines meet in the temperature-doping phase diagram. Neglecting the fluctuations of 
the magnetic field and the quenched randomness introduced by doping, see, e.g., Ref. [19| 



for a critical discussion of this point, one may consider the 0(3)©0(2)-symmetric LGW 
Hamiltonian to infer the critical behavior at the MCP, see, e.g., Refs. |8^j3^JT^j26|JT^JT5|1 . In 



particular, the analysis of Ref. ]3"3{] , which uses the projected SO (5) model |ST| as a starting 
point, shows that one can use Eq. ( |1.1| ) as an effective Hamiltonian. Different scenarios 
have been proposed for the critical behavior at the MCP. In Refs. ||[1^,|14|], it was specu- 
lated that the MCP is a bicritical point where the 0(5) symmetry is asymptotically realized. 
On the other hand, on the basis of the 0(e) results of Refs. Pj3j, Refs. [0)111 predicted a 
tetracritical behavior governed by the BFP. However, since it was expected that the BFP is 
close to the 0(5) FP, it was suggested that at the MCP the critical exponents were in any 
case close to the 0(5) ones. 

The 0(5)-symmetric scenario would require the stability of the 0(5) FP. Evidence in 
favor of this picture has been recently claimed using Monte Carlo simulations for a five- 
component 0(3)©0(2)-symmetric spin model [^,20]. The numerical results show that, 
within the parameter ranges considered, the scaling behavior at the MCP is consistent with 
an 0(5)-symmetric critical behavior. Similar results have been obtained in Ref. [[35] by a 
quantum Monte Carlo study of the quantum projected SO (5) model in three dimensions. On 
the other hand, the interpretation of these numerical results as an evidence for the stability of 
the 0(5) FP fll4l , |20| is untenable, because the results discussed in this paper definitely show 
that the 0(5) FP is unstable, and that the asymptotic approach to the MCP is characterized 
by a decoupled critical behavior or by a first-order transition. The 0(5) symmetry can be 
asymptotically realized only by tuning a further relevant parameter, beside the double tuning 
required to approach the MCP. We note that the crossover exponent <j)^4, related to the spin- 
4 pertubation of the 0(5) FP, « 0.15, is much larger than its 0(e) approximation, i.e., 
04,4 ~ ^re from which one would obtain 04^ »s 0.04 setting e = 1. It is of the same order of 
the crossover exponent appearing in many other physical systems. For instance, in randomly 
dilute uniaxial magnetic materials — a class of systems whose asymptotic critical behavior 
has been precisely observed both numerically and experimentally, see, e.g, Refs. [44 , ^2] for 
reviews — the pure Ising fixed point is unstable with crossover exponent <fi ~ 0.11, which is 
even smaller than the above-reported estimate for the 0(5) case. Therefore, contrary to 
some recent claims ||35|| , it cannot be excluded a priori that experiments are able to observe 
the unstable flow out of the 0(5) fixed point, even in those systems with a moderately small 
breaking of the 0(5) symmetry, such as the projected SO (5) model. Evidence in favor of 
a tetracritical behavior has been recently provided by a number of experiments, see, e.g., 
Refs. p7|j28|j29|j30|j3T|j3^] , which seem to show the existence of a coexistence region of the 
antiferromagnetic and superconductivity phases. The possible coexistence of the two phases 
has been discussed in Refs. f?[, fL5| , jS3|l . 

Finally, we mention that a multicritical behavior with two XY order parameters is ex- 
pected in liquid crystals, at the nematic-smectic-A-smectic-C multicritical point |T7j and 
in the presence of ferromagnetic and nematic interactions []34 l. 
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